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Abstract 



^ ■ Because of the limited size of the satelHte-borne instruments, it has not been 

On , possible to observe the flux of gamma ray bursts (GRB) beyond GeV energy. 

f— ^ \ We here show that it is possible to detect the GRB radiation of TeV energy 

CN \ and above, by detecting the muon secondaries produced when the gamma rays 

^ST ' shower in the Earth's atmosphere. Observation is made possible by the recent 

0^ ■ commissioning of underground detectors (AMANDA, the Lake Baikal detector 
and MILAGRO) which combine a low muon threshold of a few hundred GeV 

*pH| or less, with a large effective area of 10^ m^ or more. Observations will not 

O ' only provide new insights in the origin and characteristics of GRB, they also 

-*-? ' provide quantitative information on the diffuse infrared background. 

> 
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I. MUON BURST ASTRONOMY. 

High energy gamma rays produce muons when interacting in the Earth's atmosphere. 
These can be efficiently detected, and the direction of the parent gamma ray reconstructed, in 
relatively shallow underground "neutrino" detectors [||] such as the now operating Antarctic 
Muon And Neutrino Detector Array (AMANDA) and Lake Baikal telescope [^. These 
instruments are positioned at a modest depth of order 1 kilometer and are therefore sensitive 
to muons with energies of a few hundred GeV, well below the TeV thresholds of other deep 
underground detectors such as Superkamiokande and MACRO [|l|. They are therefore able to 
detect muons from primary gamma rays of TeV energy and above, with a very large effective 
telescope area of 10^ m^, or more. They infer the photon direction by reconstructing the 
secondary muon track with degree-accuracy. Although muons produced by gamma rays from 
astronomical sources compete with a large background of atmospheric cosmic ray-induced 
muons, during the short duration of a GRB this background is manageable. Using the time 
stamp provided by satellite observation, the detector only integrates background over the 
very short time of the burst which is of order one second. Unlike air Cherenkov telescopes, 
muon detectors cover a large fraction of the sky with a large duty cycle, e.g. essentially 100% 
efficiency for more than a quarter of the sky in the case of the AMANDA detector with a 
South Pole location. 

In this paper we demonstrate how large area detectors operating with a few hundred GeV 
muon threshold, or less, provide a unique window of opportunity for observing GRB. While 
the fluxes of TeV photons are reduced by one or more orders of magnitude compared to GeV 
photons observed with satellites, and while only one percent of the gamma rays will produce 
a detected secondary muon, observation of GRB is possible because the detectors are four 
orders of magnitude larger than, for example, the EGRET instrument on the Compton 
Gamma ray Observatory |^ . 

By the most conservative estimates, we predict order 1 muon per year correlated in 
time and direction with GRB for the operating detectors, and as high as 50 per year for 
the most conservative GRB flux estimate normalized to the observed diffuse GeV cosmic 
background [QQ . The unknown energetics of GRB above GeV may yield much higher rates; 
see, for instance, Ref. Pj. Interestingly, the muon count is similar for a single nearby burst 
at redshift z = 0.1, and a flux of a few muons or more in a 1 second interval in coincidence 
with a gamma ray burst cannot be missed. For a cosmological distributions such an event 
occurs within 2 ~ 3 years of observation, taking into account all observational constraints. 

Failure to observe a signal within a few years would establish a cutoff on the GRB flux not 
much above the GeV- sensitivity of the operating detectors, or point to an unexpectedly large 
infrared diffuse background absorbing the TeV gamma rays over cosmological distances. The 
two possibilities can be distinguished on the basis of observation of the occasional nearby 
burst. There is also the possibility of cosmological evolution of the sources. The observations 
are obviously relevant to the proposal that GRB are the sources of the highest energy cosmic 
rays [§. 

Also of interest here is the MILAGRO detector [Q and the proposed HANUL experiment 
in Korea P]. Although of more modest size compared to "neutrino" telescopes, MILAGRO's 
muon threshold is only 1.5 GeV because of its location at the surface. Sensitivity to gamma 
rays of lower energy results in an increased flux, and compensates for its smaller effective 



area. We should here point out that MILAGRO, as well as the other instruments discussed, 
have other capabilities to study GRB. For instance, AMANDA has sensitivity to the MeV 
neutrinos produced in the initial collapse |T^, and to neutrinos of TeV energy, and above 
0. The MILAGRO detector, beyond counting muons, can efficiently reconstruct gamma 
ray showers once their energy exceeds hundreds of GeV M . 



II. MUONS IN GAMMA RAY SHOWERS 



When gamma rays interact with the atmosphere they initiate cascades of electrons and 
photons, but also some muons. The dominant source of muons is the decay of charged pions 
which are photoproduced by shower photons. The number of muons with energy above i?^ 
in a shower initiated by a photon of energy E^ was computed some time ago ||1T|. For E^j^ 
in the range 0.1 to 1 TeV the number of muons in a photon shower can be parametrized as 



N^{E^, > E, 



2.14 X 10-^ 



cos 6 [E^l cos 6) 



E^ 



{E^/cos( 



(1) 



with energy in TeV units. The parametrization is valid for E^/ E^ > 10. This parametriza- 
tion is adequate to calculate muon rates in the AMANDA and Baikal detectors since both 
have muon thresholds of the order of a few hundreds GeV. 

For muons energies below 0.1 TeV, muon decay and muon energy losses in the atmosphere 
must be taken into account. To accomplish this we have used a linear shower Monte Carlo 
simulation which follows all shower particles down to muon threshold. It accounts for all 
decay modes and energy losses ||12[. We have calculated the number of muons at sea level 
with energy above 1.5 GeV, which is MILAGRO's threshold for muon detection. We have 
parametrized the results as a function of initial photon energy at different zenith angles for 
use in subsequent calculations. 

The Monte Carlo calculates the production of muons by the decay of photoproduced 
pions. The photoproduction cross section is obtained by interpolation of 7-proton data, 
including the most recent high energy measurements [|T^] performed with accelerators. It 
is converted to a 7-air cross section using the empirical A^'^^ dependence on the atomic 
number. A constant diffiactive cross section of 0.194 mb proceeding via p production has 
been included in the Monte Carlo simulation. For high energy muons, the Monte Carlo 
results are adequately parametrized by Eq. (1). 

We here neglected the production of muons by direct /i-pair production and leptonic 
decay of charmed particles which only contribute to the flux of muons with energy above 



several TeV HIT 



III. GRB MUON RATES FROM ENERGY CONSIDERATIONS 



It is thought that GRB are produced when a highly relativistic shock with Lorentz factor 
of order 100 or higher dissipates its kinetic energy in collisions with the inter-stellar medium 
or in internal collisions within the relativistic ejecta. The observed 7-rays are most likely 
produced by synchrotron radiation of electrons accelerated in the shock, possibly followed 



by inverse Compton scattering, when the relativistic shell becomes optically thin to pair 
production |]14| . 

The observed injection rate in the Universe is -Egrb = 4 x 10*^^ erg Mpc~'^ yr"^- Inter- 
estingly, this is equal to the injection in cosmic rays beyond the ankle in the spectrum near 
10® TeV [0. The BATSE satellite observes, on average, 1 GRB per day with an efficiency 
of about 25%. Therefore the average energy emitted per burst, -Egrb, is : 

'^i ( D V ( E \ /4day"^\ 

i?OKB - 3 X 10 ergs (^^53^^^ j ( 4 x lO^^rg Mpc-3 yr-i j (^^ j " ^'^ 

R is the burst rate. Throughout our analysis we will assume a cosmological distribution of 
GRB, with a distance oi D = 3000 Mpc to the average burst. Assuming no beaming, this 
energy corresponds to an average fluency per GRB, Fgrb, close to the one observed: 



. ..-St -2^ 3000Mpc y^ Egrb ^ .„. 

^GRB = 3xlO Jm [—^^) ,,1052 ...J - (3) 



Although in line with observations, this may be an underestimate because present detectors 
provide no information on the energetics of GRB above GeV-energy. Some authors [^,|^ 
have, for instance, raised the interesting possibility that GRB are the origin of the diffuse 
extragalactic gamma ray background. This association requires a photon energy of order 
-Eqrb ~ 10^^-10^® ergs per burst, depending on the GRB occurrence rate assumed. Most of 
this energy is concentrated in the high energy tail of the photon spectrum, above E^ > 1 TeV, 
and is therefore not included in the energy balance of equation ^ where only observed photons 
of GeV-energy and below have been considered. There are models that can accommodate 
the large total energy required by this scenario, for instance those where a complete neutron 
star is converted into photons |l^ . 



A typical GRB spectrum exhibits a high energy tail which, above a few hundreds keV, 
can be parametrized by a power law: 

^ = — ^ 10~^2 cm"2 s"^ (A) 

dE^ E(7+i) ^^ cm s , i4j 

where energies are in TeV. Typical values for the observed spectral index 7 range from 0.8 



to 1 1 14]. The photon spectrum of the average burst is normalized by energy conservation: 



^^"'' dE E ^= ^^^^ , (5) 

where At is the average duration of the burst of order 1 second. -E'^min {E^ma.x) is the 
minimum (maximum) energy of the photons emitted by the burst. For -E^min we will use 
1 MeV throughout. Observations indicate indeed that a negligible fraction of the energy is 
emitted in X-rays. 

The TeV flux obtained by extrapolation, is too low to be detected by the present satellite 
experiments because of their small telescope area. The high energy behavior of the GRB 
spectrum is therefore a matter of intense speculation. Bursts where isolated photons reach 
tens of GeV energy have been detected, suggesting the extension of the spectrum to the TeV 



range [14|. As previously discussed, the extension of the photon spectrum to TeV energy 



with a rather flat spectral index (7 ~ 0.5) is an essential feature of models accommodating 



the diffuse GeV background. Also, the HEGRA group [|I^ has observed an excess of gamma 
rays with E^ >16 TeV in temporal and directional coincidence with GRB 920925c. 
The value of F^ obtained through above normalization procedure is given by: 



F, = 3xlO'''™^'P'=^ ' ^™= 




3 X 10^^ ergs 



- <~;n"^]-^ for 7 7^ 1 
for 7 = 1 



(6) 



Note that, to a first approximation, the normalization factor F^ only depends on £"^111111 
(-E^max) for 7 > 1 (7 < 1). The reason for this is clear, when 7 > 1 the spectrum is very 
steep and most of the GRB energy is concentrated in lower energy photons. The situation 
is reversed in the case of 7 < 1. We will explore the dependence of the muon rate on 7 
further on. 

Absorption of gamma rays in the infrared, optical and microwave backgrounds is a de- 
termining factor in the gamma ray flux observed at Earth. The mean free path of a TeV 
photon in the diffuse background is thought to be less than a few hundred Mpc, although 
this estimate is very uncertain because of our poor knowledge of the diffuse infrared back- 
ground. Nevertheless, for a cosmological population of astrophysical objects such as GRB, 
absorption is expected to reduce the detected flux of TeV photons. The detection could 
therefore be dominated by the closest bursts. 

We first compute N^{E^,> E^), the number of muons with energy in excess of E^, 
produced in a photon shower of energy E^. The muon rate is subsequently derived by 
convolution with the gamma ray spectrum, 

N,{> E,) = dE, N,{E„ >E,)—^. (7) 

Here -E^th? the minimum photon energy needed to produce muons of energy E^, is given by 
Eryt^ ~ 10 X E^/cos6, where 6 is the zenith angle at which the source is observed. 

Our results imply that the operating AMANDA, Baikal and MILAGRO detectors can 
do GRB science in an energy range not covered by astronomical telescopes. While detection 
may be marginal for the most conservative estimates of the high energy GRB flux, this will 
not be the case for future detectors on the drawing board such as AMANDA II, IceCube 



and Antares [17 



A. AMANDA and Lake Baikal as Gamma Ray Telescopes 

The vertical threshold for muon detection in the AMANDA detector at a 1.5 km depth 
is ~ 350 GeV. For the shallower depth of lake Baikal this threshold is ~ 150 GeV. For 
threshold energies in this range we may make use of parametrization (|l]) which is adequate 
for E^ between 0.1 and 1 TeV [|TB|. As previously mentioned, absorption of gamma rays in 



cosmological backgrounds has to be taken into account, as well as the source distribution of 
GRB in redshift z: 



N,{> E,){jT-') = A{m') At-t / / dz dE.—L, e'^^^-^) 

i JO JO JE^n^in Ey 

where i?^ is the vertical threshold energy of the detector, in TeV units. r^E^, z) is the optical 
depth of a photon with energy E^ originating at a distance z. For illustration, we will use 



the diffuse photon background of reference |T^ throughout. dRo^isldz is the cosmological 



distribution of GRB, i.e. the number of GRB per unit time and distance z. Following pO|JT4 
we will assume cosmologically distributed standard candles with no source evolution. The 
GRB distribution has been normalized to 365 events per year observed by BATSE below a 



redshift of Zmax = 2.1 ||2^. A is the effective area of the detector. The dependence of F^ on 
the characteristics of the burst is given by Eq.(|^), which can be used for scaling purposes. 

AMANDA and Baikal observations require a time-stamp for background rejection, they 
can therefore only see bursts previously observed by satellite experiments. The BATSE rate 
of 1 burst/day corresponds to an isotropic distribution on the sky of 365 bursts/year. For the 
most conservative estimate with -E^max = 10 TeV and a "standard" burst with total energy 
-^GRB = 3 X 10^^ ergs and a spectral index 7 = 1 , we predict 1 muon per year correlated 
in time and direction with GRB for the AMANDA telescope, and 89 per year for the most 
conservative GRB flux estimate normalized to the observed diffuse GeV cosmic background 
PJ^. The unknown energetics of GRB above GeV may yield much higher rates; see, for 
instance ref H. The results are subject to the uncertainty associated with the absorption 



on infrared light which we computed following Stecker |jT9[. The rates are essentially the 
same for the Lake Baikal experiment. 

Interestingly, the muon counts are similar for a single nearby burst at redshift 2; = 0.1, 
and a flux of a few muons or more in a 1 second interval in coincidence with a gamma 
ray burst provides a striking experimental signature which is hard to miss. In Figs. 1 and 
2 we plot the number of muons per burst versus spectral index 7 for a burst at 2; = 0.1. 
We assume again that the photon spectrum only extends to -E^max = 10 TeV with a total 
energy £'grb = 3 x 10^^ ergs. We assumed an effective area oi A = 10^ m^ for AMANDA 
and 10^ m^ for the Lake Baikal detector [0,0. The results are shown with and without 
absorption on infrared light. Also the dependence of the rate on zenith angle 9 is illustrated. 
The rapid decrease of the number of events when 9 increases is a consequence of the increase 
of the energy threshold because muons have to penetrate an increasing amount of matter to 
reach the detector. The smaller area of Baikal is compensated by its lower muon threshold, 
especially for larger zenith angles. 

In Fig. 3 we show the rate per burst in AMANDA detector as a function of redshift for 
a model with spectral index 7 = 1 and total energy -Egrb = 3 x 10^^ ergs. Also shown is 
the flux for energetics which accommodates the diffuse GeV background. Absorption has 
been included, as in Figs. 1 and 2. Shown as a dashed line is the number of days towards a 



burst at the corresponding redshift assuming the GRB distribution of reference ||2^ 



An important question is whether the signal is observable given the large background of 
muons of cosmic ray origin penetrating the detectors. In this respect the great advantage 
of GRB detection is that the detector only integrates background for the short duration of 
the burst, typically 1 second. Furthermore, the background can be limited to a circle in the 
sky of angle 59 around the direction of the burst. Here 59 is the angular resolution of the 

6 



detector which is typically a few degrees. It can be sharpened up by quality cuts, but this 
inevitably results in a loss of effective area. The muon background intensity from cosmic 
ray showers I^{0) has been measured [0,0. The muon background at zenith 6 is given by: 



Noise = I^{9)x Ax 69^ 
The signal-to-square root of noise then scales as 



N 



69 



(9) 



(10) 



which simply expresses that the sensitivity is improved for larger area and for better angular 
resolution. Note that our search is not sensitive to photons produced over time-scales much 
larger than seconds, for instance to photons possibly produced by propagation of the external 
shock in the interstellar medium. 

In Table I we present the number of muons per burst detected by AMANDA. The first 
column is the redshift of the burst, the second (third) column is the signal without (with) 
absorption of gamma rays in intergalactic backgrounds. In each columns the results is given 
for zenith angle ^ = 60 and ^ = 0. The fourth column gives the signal to square root of 
noise corresponding to column 3. The 5th column shows the time in days between bursts 
with z lower than the value shown in the first column. We here assumed a "standard" burst 



with spectral index 7 = 0.8 and E, 



7max 



10 TeV. 



0.05 
0.1 
0.5 
1.0 



Muons 
per burst 



5-25 

1 -6 
0.05 - 0.25 
0.01-0.05 



Muons 
per burst 



2-12 

0.06-0.6 

5 X 10-1° - 3 X 10-6 

- 3 X 10-11 



Signal to 

Vn 



oise 



60 133 

se 59 

1.8 6.6 

se 59 

1.5x10-** 3.3x10-5 







3.3x10-™ 



Time 

(days) 



8400 

1200 

24 



Table I: Event rates in AMANDA detector. For models accommodating the diffuse 
GeV background the rates are larger by 2 orders of magnitude. The first column 
is the redshift of the burst, the second (third) column is the signal without (with) 
absorption of gamma rays in intergalactic backgrounds. In each columns the results 
is given for zenith angle 6' = 60 and 9 = 0. The fourth column gives the signal to 
square root of noise corresponding to column 3. The 5th column shows the time in 
days between bursts with z lower than the value shown in the first column. 

The role of absorption for gamma ray detection is evident. Nearby bursts with z < 0.1 
provide the best opportunity for detection. Their frequency is of the order of 1 burst every 
2 ~ 3 years. 



B. The MILAGRO Telescope 



The muon energy threshold for MILAGRO detector is only 1.5 GeV. It has an effective 
area A = 1.5 x 10^ m^ and an intrinsic angular resolution of about 3°. This requires the use of 



a 4.7° bin-size to collect ~ 70% of the gamma ray events in which the muon background from 
cosmic rays is ~ 900 Hz [§ . In Table II we present the number of muons per burst observed 
by MILAGRO for two maximum energies and two spectral slopes of the photon spectrum. 
The results are shown corresponding to the "standard" burst energy, and to the energy 
required in the models where GRBs are responsible for the diffuse extragalactic gamma ray 
background. The absorption of the photon flux has not been taken into account. MILAGRO 
with a threshold for muon detection as low as 1.5 GeV is sensitive to photons with energy 
~15 GeV which are not absorbed in the intergalactic backgrounds. Therefore detection of 
bursts at large redshifts is possible. For a spectral index 7 = 1, equal amounts of energy 
are stored in every logarithmic interval. The number of low energy photons will be larger 
than high energy ones by the ratio of energies. On the other hand, the number of muons 
in a photon shower scales roughly with its energy. This will compensate the contribution 
to the number of muons from low and high energy photons if it were not for absorption 
which predominantly reduces the high energy part of the photon spectrum. Counting of low 
energy muons may therefore be more efficient than reconstruction of the photon shower, the 
conventional MILAGRO method, because it raises the threshold to hundreds of GeV. This 
is especially true for the detection of the average {z ~ 1) GRB. 
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Spectral Slope 


-t^7max 


Signal 


Signal to 


(ergs) 


7 


(TeV) 


per burst 


VNoise 


3 X 10^2 


1.0 


10. 


24 


0.82 


3 X 10^2 


0.8 


10. 


57 


1.9 


3 X 10^2 


0.5 


10. 


95 


3.2 


3 X 10^2 


1.0 


0.1 


2.3 


0.08 


3 X 10^2 


0.8 


0.1 


5.0 


0.16 


3 X 10^2 


0.5 


0.1 


9.2 


0.3 


1054 


1.0 


10. 


785 


26.4 


1054 


0.8 


10. 


1920 


64.0 


1054 


0.5 


10. 


3230 


107.4 


1054 


1.0 


0.1 


75 


2.5 


1054 


0.8 


0.1 


165 


5.6 


1054 


0.5 


0.1 


305 


10.0 



Table II: Sensitivity of MILAGRO to a GRB at zenith angle 6=0, for several assump- 
tions on total energy, spectral index and maximum energy. Absorption of gamma rays 
is not included. 
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FIG. 1. Number of muons in AMANDA per "standard" burst (see text) at 2; = 0.1 as a function 
of the spectral index of GRB photon spectrum. Curves are shown for different zenith angles of 
observation and with and without absorption of photons in the infrared background. 
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Figure 2 
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FIG. 2. Same as Fig. 1 for Lake Baikal detector. 
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Figure 3 
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FIG. 3. Number of muons per burst in AMANDA as a function of redshift. Results include 
absorption by the infrared background and are shown for a "standard" burst and a burst normalized 
to the observed GeV diffuse gamma ray background. The dashed line shows the time in days until 
a burst occurs in the field of view of AMANDA, at a z value lower than the one indicated on the 
horizontal axis. 



12 



